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^ ! Abstract 

oo : 

^ \ We explore 'benchmark surfaces' suitable for studying the phenomenology of Higgs bosons 

^ , in the minimal supersymmetric extension of the Standard Model (MSSM), which are chosen 
^ ■ so that the supersymmetric relic density is generally compatible with the range of cold dark 
^ ■ matter density preferred by WMAP and other observations. These benchmark surfaces are 
^ \ specified assuming that gaugino masses mi/2, soft trilinear supersymmetry-breaking param- 
^ , eters Aq and the soft supersymmetry-breaking contributions mo to the squark and slepton 
■ masses are universal, but not those associated with the Higgs multiplets (the NUHM frame- 
work). The benchmark surfaces may be presented as (M^, tan (5) planes with fixed or system- 
atically varying values of the other NUHM parameters, such as mo, mi/2, Aq and the Higgs 
mixing parameter /z. We discuss the prospects for probing experimentally these benchmark 
surfaces at the Tevatron collider, the LHC, the ILC, in B physics and in direct dark-matter 
detection experiments. An Appendix documents developments in the FeynHiggs code that 
enable the user to explore for her/himself the WMAP-compliant benchmark surfaces. 
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1 Introduction 



Some of the best prospects for probing the minimal super symmetric extension of the Standard 
Model (MSSM) [1,2] might be offered by searches for the bosons appearing in its extended 
Higgs sector. It may be challenging to distinguish between the lightest MSSM Higgs boson 
and a Standard Model (SM) Higgs boson with the same mass, and searches for MSSM 
Higgs bosons are, in many ways, complementary to searches for supersymmetric particles as 
avenues to establish the existence of physics beyond the SM. 

Searches at the Tevatron collider are closing in on the possible existence of an SM- 
like Higgs boson over a limited range of low masses [3-5] , and are also starting to encroach 
significantly on the options for heavier MSSM Higgs bosons, particularly at large tan f3 [6-10]. 
Studies have shown that experiments at the LHC will be able to establish the existence or 
otherwise of an SM-like Higgs boson over all its possible mass range, and also explore many 
options for the heavier MSSM Higgs bosons [11-14]. On the other hand, the LHC might well 
be unable to distinguish between the lightest MSSM Higgs boson and an SM Higgs boson of 
the same mass. The ILC would have better chances of making such a distinction [15-22], and 
might also be able to produce the other MSSM Higgs bosons if they are not too heavy [15-19]. 
CLIC would also be able to study a light SM-likc Higgs boson, as well as extend the search for 
MSSM Higgs bosons to much higher masses [23]. Searches for new phenomena in B physics, 
including rare decays such as 6 — > sj, Bg fJ^'^IJ'~ and Bu —>■ ru, also have good potential to 
explore the MSSM Higgs sector and, at least in some specific MSSM scenarios, electroweak 
precision observables (EWPO) may also provide interesting constraints [24,25]. In parallel to 
these accelerator searches for MSSM Higgs bosons and their effects, non-accelerator searches 
for supersymmetric dark matter [26,27] will also be able to explore significant regions of the 
MSSM Higgs parameter space [28-30], since the exchanges of massive MSSM Higgs bosons 
have significant impacts on dark matter scattering cross sections. 

In order to correlate the implications of searches at hadron colliders and linear colliders, 
in B physics, in dark matter searches and elsewhere, it is desirable to define MSSM Higgs 
benchmark scenarios that arc suitable for comparing and assessing the relative scopes of 
different search strategies, see, e.g., Refs. [31-38]. 

Since the MSSM Higgs sector is governed by the two parameters (or Mfj±) and tan/9 
at lowest order, aspects of MSSM Higgs-boson phenomenology such as current exclusion 
bounds and the sensitivities of future searches are usually displayed in terms of these two 
parameters. The other MSSM parameters enter via higher-order corrections, and are conven- 
tionally fixed according to certain benchmark definitions [31-34]. The benchmark scenarios 
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commonly used in the literature encompass a range of different possibilities for the amount 
of mixing between the scalar top quarks, which have significant implications for MSSM Higgs 
phenomenology, and also include the possibihty of radiatively-induced CV violation. The 
best- known example is the so-called "m™^ scenario" [31-33], which allows the search for 
the fight CP-even Higgs boson to be translated into conservative bounds on tan /3 for fixed 
values of the top-quark mass and the scale of the supersymmetric particles [39] . The existing 
benchmark scenarios designed for the MSSM Higgs sector are formulated entirely in terms 
of low-scale parameters, i.e., they are not related to any particular SUSY-breaking scheme 
and make no provision for a possible unification of the SUSY-breaking parameters at some 
high mass scale, as occurs in generic supergravity and string scenarios. 

In applications of the existing benchmark scenarios for the MSSM Higgs sector [31-34], 
one is normally concerned only with the phenomenology of the Higgs sector itself. Besides 
the direct searches for supersymmetric particles, other constraints arising from EWPO, B- 
physics observables (BPO) and the possible supersymmetric origin of the astrophysical cold 
dark matter (CDM) are not usually taken into account. This may be motivated by the 
fact that the additional constraints from EWPO, BPO and CDM can depend sensitively 
on soft-supersymmetry breaking parameters that otherwise have minor impacts on Higgs 
phenomenology. For example, the presence of small fiavour-mixing terms in the MSSM La- 
grangian would severely affect the predictions for the BPO while leaving Higgs phenomenol- 
ogy essentially unchanged (see also Ref . [36] for a discussion of this issue) . 

In this paper we follow a different approach and adopt specific universality assumptions 
about the soft SUSY-breaking parameters, restricting our analysis of the MSSM to a well- 
motivated subspace of manageable dimensionality. It is frequently assumed that the scalar 
masses mo are universal at some high unification scale, as are the gaugino masses mi/2 
and the trifinear parameters Aq, a framework known as the constrained MSSM (CMSSM). 
In such a scenario, the heavier MSSM Higgs boson masses are fixed in terms of the input 
parameters and tan /3, so that Ma is not an independent parameter, and consequently this 
scenario is too restrictive for our purposes. However, there is no good phenomenological or 
theoretical reason why the soft supersymmetry-breaking contributions to the Higgs masses 
should not be non- universal, a scenario termed the NUHM [40-42]. Within the NUHM, 
Ma and // can be treated as free parameters for any specified values of mo,mi/2, Aq and 
tan /3, so that this scenario provides a suitable framework for studying the phenomenology 
of the MSSM Higgs sector. Since the low-scale parameters in this scenario are derived from 
a small set of input quantities in a meaningful way, it is of interest to take into account other 
experimental constraints. 
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The main purpose of this paper is to explore new benchmark surfaces for MSSM Higgs 
phenomenology that are compatible with the cosmological density of cold dark matter in- 
ferred from a combination of WMAP and other observations [43]. While in the CMSSM 
only narrow strips in (mi/2,'n^o) planes are compatible with WMAP et al. [44,45] for given 
values of Aq and tan (3, the NUHM offers the attractive possibility to specify {Ma, tan (3) 
planes such that essentially the whole plane is allowed by the constraints from WMAP and 
other observations [25]. This is done assuming that R parity is conserved, that the lightest 
supersymmetric particle (LSP) is the lightest neutralino Xi? and that it furnishes most of 
the cold dark matter required [46]. As we discuss in more detail below, compatibility with 
WMAP et al. cannot be maintained while keeping all the other NUHM parameters fixed. 
Accordingly, we discuss two examples of WMAP-compliant benchmark surfaces that are 
specified for fixed mo, /i and Aq = but varying mi/2, and two surfaces that are specified 
for fixed mi/2, rno and Aq = but varying fi. For the first two benchmark surfaces, a simple 
linear relation between mi/2 and Ma is imposed as the {Ma, tan (3) plane is scanned, whereas 
for the other two surfaces /i is varied through a relatively narrow range. 

Following the specifications of these NUHM benchmark surfaces, we then explore the pos- 
sibilities for studies of the MSSM Higgs bosons and other supersymmetric signatures across 
these {Ma, tan P) planes. We consider the electroweak precision observables, principally 

= j{g — 2)fj_ and Mh, prospects for the search for H/A t^t^ at the Tevatron, prospects 
at the LHC - including searches for /i ^ 77 and r~^r~, H/A t^t^ and r^z/, and 

measurements of the ratio oi h ^ t^t^ and WW* branching ratios, prospects at the ILC - 
including ways of distinguishing between the light MSSM h boson and an SM Higgs boson of 
the same mass by measuring (ratios of) branching ratios, prospects in B physics - including 
Bg fJ'^fJ' , b s'-f and Bu —>■ tu, and the direct detection of supersymmetric cold dark 
matter. In an Appendix we introduce developments in the FeynHiggs code that enable the 
user to explore for her/himself the WMAP-compliant benchmark surfaces. These include 
the concept of a FeynHiggs record, a new data type that captures the entire content of a 
parameter file in the native format of FeynHiggs. 

2 Specification of the Benchmark Surfaces 

As an introduction to the specification of the benchmark surfaces in the NUHM, we first 
consider a generic {Ma, tan P) plane for fixed mi/2,mo, Aq and fi, adapted from Ref. [47]. 
As we see in Fig. [Tt^a), in the {Ma, tan (3) plane for mi/2 = 600 GeV, ttiq = 800 GeV, 
fi = 1000 GeV and Aq = 0, the relic LSP density satisfies the WMAP constraint only in 
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narrow, near-vertical (pale blue) shaded strips crossing the plane. These lie to either side of 
the vertical (purple) line where m^o = Ma/^. Within the narrow unshaded strip straddling 
this line, the relic density is suppressed by rapid direct-channel annihilations to a value 
below the lower limit of the range for the cold dark matter density indicated by WMAP et 
al. This strip would be acceptable for cosmology if there were some additional component of 
cold dark matter. Outside the shaded WMAP-compatible strips, at both larger and smaller 
values of Ma, the rehc LSP density is too high, and these regions are unacceptablj^]. 

It is clear from this example that one may arrange for the relic LSP density to remain 
within the preferred WMAP range over (essentially) the entire [Ma, tan (3) plane if one 
adjusts mi/2 continuously as a function of Ma so as to remain within one of the narrow 
WMAP strips as Ma increases. Accordingly, we study a benchmark {Ma, ta.n (3) plane PI 
with the same values of rriQ = 800 GeV, fi = 1000 GeV and Aq = 0, but with varying mi/2 ~ 
|M^. Since we evaluate observables using a discrete sampling of the NUHM parameter 
space, we consider values of mi/2 lying within the small range of this central value: 

^Ma - 12.5 GeV < mi/2 < ^Ma + 37.5 GeV. (1) 
8 8 

The observables that we study do not vary significantly as mi/2 is varied across this range. 
Specifically, we use the mi/2 that gives the value of the cold dark matter density that is 
closest to the central value within the allowed range, 0.0882 < ^cdm^^ < 0.1204 [43] (see 
below) . 

Previous analyses of the CMSSM indicated that values of mi/2 and mo below 1 TeV are 
preferred, in particular by the EWPO [25,48,49] (see also Ref. [50]). Accordingly, we study 
also a benchmark {Ma, ta.n j3) plane P2 with the fixed values mo = 300 GeV, fi = 800 GeV 
and Aq = 0, with mi/2 ~ 1.2M^ again varying continuously across the plane so as to 
maintain the WMAP relationship with Ma- As before, because of our discrete sampling of 
the NUHM parameter space, we consider values of mi/2 lying within a small range of this 
central value: 

1.2Ma - 40 GeV < mi/2 < 1.2Ma + 40 GeV. (2) 

Again, the observables that we study do not vary significantly as mi/2 is varied across this 
range. 

More examples could be chosen with different fixed values of mo, fi and Aq but, as long 
as mi/2 is the parameter being varied to keep the LSP density within the WMAP range. 



^We note in passing that the LEP lower hmit on Mh excludes a strip of this plane at low Ma and/or 
tan/3 indicated by the dash-dotted (red) hne, that (pink shading) prefers relatively large tan/3 > 36, that 
6 ^ S7 excludes a (green shaded) region at low Ma and tan /3, and that the other BPO disfavour a region 
at low Ma and high tan/3 (not shown). 
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Figure 1: Sample NUHM parameter planes with two parameters varied and the other four 
fixed, adapted from Refs. [47,51]. The left plot displays a {Ma, tan (3) plane with mi/2 = 
600 GeV, mo = 800 GeV, /i = 1000 GeV and Aq = 0. The range of cold dark matter density 
preferred by WMAP and other observations is attained in two narrow (pale blue) strips, one 
on either side of the vertical solid (blue) line where m^o = Ma/2. The dark (green) shaded 
region at low Ma and low tan j3 is excluded by 6 ^ 57, and the medium (pink) shaded region 
at tan/5 > 36 is favoured by a^. The region below the (red) dot-dashed fine is excluded by 
the LEP bounds on Mh. The right plot displays a (/i. Ma) plane with mi/2 = 500 GeV, 
mo = 1000 GeV, tan/5 = 35 and Aq = 0. Here the WMAP range of cold dark matter 
density is attained in two narrow strips at roughly constant positive and negative values of 
/i, which are swept apart by rapid annihilation when Ma ~ 2m^o. The dark (green) shaded 
region at ;U < is excluded by 6 ^ 57, and the < /i < 760 GeV strip (pink shading) is 
favoured by a^. The region below the (red) dot-dashed line again is excluded by the LEP 
bounds on Mh, and the region between the vertical (black) dashed lines has m-± < 104 GeV. 



a similar relationship between mi/2 and Ma will always apply. The only fiexibility in the 
choice of mi/2 is whether one wishes to stay within the left or right near-vertical shaded 
strip. However, the corresponding values of mi/2 do not differ greatly, and neither do the 
corresponding phenomenological signatures, though the lightest Higgs boson mass can be 
somewhat sensitive to this choice. The values of mo and (to a lesser extent) fi have far more 
impact on the phenomenology, and the benchmark choices we have made: mo = 800 GeV 
for PI and mo = 300 GeV for P2, provide significant and interesting differences worthy of 
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examination. 

We also study two other {Ma, tan j3) planes, whose motivation can be gained from ex- 
amination of the (/i, Myi) plane shown in Fig. [T](b), which is adapted from Ref. [51]. We 
see that, for a fixed choice of values of mi/2 = 500 GeV, rriQ = 1000 GeV and Aq = 0, 
there is a narrow strip of values of yU ~ 300 —>■ 350 GeV where the relic density lies within 
the WMAP range for almost all values of Ma- The exception is a narrow strip centred on 
Ma ~ 430 GeV, namely the rapid-annihilation funnel where m^o ~ Ma/2, which would be 
acceptable if there is some other source of cold dark matter. This funnel is narrower (wider) 
for smaller (larger) values of tan /3, but its location in /i does not vary much as a function of 

Motivated by this example, we explore two benchmark surfaces with different fixed values 
of mi/2 and tuq, and fi varying within a restricted range chosen to maintain the LSP density 
within or below the WMAP range. The first example of such a benchmark plane, P3, has 
fixed mi/2 = 500 GeV, mo = 1000 GeV and Aq = 0, with /i in the range 

= 250 - 400 GeV. (3) 

In the following, we evaluate observables for a discrete sampling of /i values within this range. 
Since the corresponding variation of the particle mass spectrum is quite small, the impact 
of the variation of /i on the observables discussed below is negligible. 

The other example of such a benchmark plane, P4, has fixed mi/2 = 300 GeV, mo = 
300 GeV and Aq = 0, with in the range 

/X = 200 - 350 GeV. (4) 

As in the previous case, the LSP density lies within the WMAP range except for a small 
range of Ma ~ 2m^o where the density is below the preferred range. However, again this is 
acceptable if there is some other component of cold dark matter. The parameter choices for 
this and the other NUHM benchmark surfaces are summarized in Tab. [^. 

A likelihood analysis of these four NUHM benchmark surfaces, including the EWPO Myy, 
sin^^eff, r^, (g - 2)^ and M^ and the BPO BR(6 ^ 57), BR(5, ^ /i+/i"), BR(5„ ^ rz/^) 
and AMb, was performed recently in Ref. [25]. The lowest value in each plane, denoted 
as Xmin5 is shown in the rightmost column of Tab. [l], corresponding to the points labeled as 

^We note in passing that the LEP lower hmit on Mh excludes a strip of this plane at low Ma indicated 
by the (red) dash-dotted line, and the LEP lower limit on the chargino mass excludes values of /i between 
the two vertical (black) dashed lines. 

A minor change in the best-fit point and the Xmin ocurred for the P2 scenario in comparison with 
Ref. [25] due to a slightly different choice of the TO1/2 values. 
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Table 1: The four NUHM benchmark surfaces are specified by the above fixed and varying 
parameters, allowing Ma and tan/5 to vary freely. All mass parameters are in GeV. The 
rightmost column shows the minimum value found in each plane at the points labelled as 
the best fits in the plots. 

the best fits in the plots below. We display in each of the following figures the locations 
of these best-fit points by a (red) cross and the Ax^ = 2.30 and 4.61 contours around the 
best- fit points in the (M^jtan/?) planes for each of these benchmark surfaces. These con- 
tours would correspond to the 68 % and 95 % C.L. contours in the (MA,tan/3) planes if 
the overall likelihood distribution, C oc e~^^/^, were Gaussian. This is clearly only approx- 
imately true, but these contours nevertheless give interesting indications on the regions in 
the (MA,tan/3) planes that are currently preferred. The varied parameter in each scenario 
(i.e. mi/2 in PI, P2 and ji in P3, P4) is chosen such that the cold dark matter density is 
closest to the central value within the allowed range, 0.0882 < Q^cnuh^ < 0.1204 [43]. 

On surfaces PI and P2, where mi/2 scales with Ma so as to remain in the funnel region, 
much of the mass spectrum scales with Ma- Specifically, the lightest neutralino and chargino 
masses simply scale in direct proportion to Ma for these surfaces. The light squark masses 
and stau masses also scale with mi/2 (and hence Ma)-, though the latter arc also slightly 
dependent on tan/? as well. In the range Ma < 1 TeV displayed in these planes, the light 
squark masses range up to ~ 2.3 TeV for surface PI, within reach of the LHC. However, 
because of the relatively large values of mg, the light squarks arc beyond the current reach 
of the Tevatron collider even at low Ma (and hence mi/2)- For P2, the light squark masses 
range up to ~ 1.7 TeV. 

Turning to surfaces P3 and P4, because they have fixed values of mi/2 and mo, there 
are very small variations in the sparticle mass spectra across these planes. For example, the 
lightest neutralino and chargino masses are determined primarily by mi/2, and so they both 
take almost constant values on the benchmark surfaces. Similarly, the light squark masses 
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are determined by a combination of mi/2 and mo and show little dependence on either Ma 
or tan (3. On the other hand, the lightest stau mass has a slight dependence on tan (3, due to 
the variable splitting of the third-generation sparticle masses. These mass splittings increase 
at large tan (3, leading to smaller stau masses. 

We display in each plane the region excluded (black shaded) at the 95 % C.L. by the 
LEP Higgs searches in the channel e"'"e~ — > Z* — Zh^ H [52,53]. For a SM-like Higgs boson 
we use a bound of > 113 GeV. The difference from the nominal LEP mass limit allows 
for the estimated theoretical uncertainty in the calculation of for specific values of the 
input MSSM parameters [54]. In the region of small Ma and large tan/3, where the coupling 
of the light CP-even Higgs boson to gauge bosons is suppressed, the bound on Mh is reduced 
to Mh > 91 GeV [52]. 

3 Electroweak precision observables 

In this Section we summarize key predictions for electroweak precision observables (EWPO) 
over the four benchmark surfaces. In Ref. [25] it was shown that Mw, sin^ 6'eff and Tz 
agree within ~ 1 cr with the current experimental value over all the benchmark surfaces. 
Since their variations are relatively small, we do not display these observables in this paper, 
though they are included in the overall function. Here we focus on two other EWPO, 
namely the mass of the lightest Higgs boson, M^, and the anomalous magnetic moment of 
the muon, = ^{g — 2)^. 

The evaluation of Mh is performed using FeynHiggs [54-57]. In Fig. [2] we show the 
contours for Mh = 113, 114, 115, 116, 117, 118 and 120 GeV. As discussed in the previous 
Section, the boundary of the region excluded by the LEP searches for the lightest MSSM 
Higgs boson does not coincide with the nominal limit Mh = 114.4 GeV on the mass of a 
Standard Model Higgs boson. Nevertheless, it can be seen in Fig. [2] that the Ax^ = 2.30 and 
4.61 contours are highly correlated with the Mh contours at low values of Ma and tan (3. This 
is a consequence of the fact that the full likelihood information from the LEP Higgs exclusion 
limit (as well as the theoretical uncertainty) is incorporated into the overall function (see 
Ref. [25]). Note that for the plane P4 (and to a lesser extent P3) the maximum value for 
the Higgs mass is limited by the relatively low value of mi/2- 

Concerning a^, we recall that, according to a recent evaluation of the Standard Model 
contribution based on low-energy e~^e~ data, there is a discrepancy with the experimental 
measurement by the E821 Collaboration [62,63]. It would be premature to regard this 
deviation as solid evidence for new physics. However, within the SUSY framework we explore 
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Figure 2: The (M^,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) P2, (c) 
P3 and (d) P4, displaying the contours of Ax^ found in a recent global fit to EWPO 
and BPO [25]. All surfaces have Aq = 0. We also display individually the contours of 
Mfi found using FeynHiggs [54-57] and the contours of found using Refs. [58-61]. The 
1(2)-(T range for is demarcated by dashed (solid) lines. The dark shaded (black) region 
corresponds to the parameter region that is excluded by the LEP Higgs searches in the 
channel e+e' ^ Z* ^ Zh,H [52,53]. 
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here, this discrepancy does impose a significant constraint on the parameter space, and makes 
an important contribution to the global function whose contours are shown in Fig. O Our 
evaluation of is based on Refs. [58-61], which yields [64,65]: 

^exp _ ^theo _ (27.5 ± 8.4) X 10-^°, (5) 

equivalent to a 3.3-cr effeciS. In Fig. [2] we show the contours Aa^ = 10.7,19.1,35.9,44.3 x 
10~^° for the net supersymmetric contribution to a^. 

In the case of surface PI, we see that the best-fit point corresponds to ~ 118 GeV 
and Attfj, ~ 10.7 x 10~^°. In most of the displayed region of the surface that is favoured at 
the global Ax^ < 4.61 level, Aa^ is considerably lower than the range favoured in eq. ([5]). 
In the case of surface P2, the best-fit point has ~ 118 GeV, and Aa^ is within the 1-a 
range given by eq. ([5]). In the case of surface P3, the best-fit point has Mh > 118 GeV and 
again a low value of Aa^. Finally, the best-fit point in surface P4 has Mh ~ 115 GeV and an 
excellent value of Aa^, according to eq. ([5]). The fact that the best-fit points do not always 
have favoured values of Aa^ refiects the importance of other precision observables, notably 
the BPO discussed later. 



4 Tevatron Phenomenology 

We first consider how experiments at the Tevatron collider in the next years could probe the 
benchmark surfaces PI, P2, P3 and P4. We consider one possible Tevatron signature for 
the MSSM Higgs sector, namely H/A t^t^, for which expectations are evaluated using 
the results from Ref. [70]. They are based on the expectation of a 30% improvement in the 
sensitivity with respect to Ref. [6]. We see in Fig. [3] that, at the Tevatron with 2 (4, 8) fb^^ 
of integrated and analyzed luminosity per experiment^, the channel H/A t^t^ would 
provide a 95% C.L. exclusion sensitivity to tan [3 ~ 35(30, 25) when Ma ~ 200 GeV, and the 
sensitivity decreases slowly (rapidly) at smaller (larger) M^. In the case of the benchmark 
surface PI, 8 fb~^ would start accessing the region with Ax^ < 4.61. For P2, however, the 
area accessible to the Tevatron is not visible in the figure since it is completely covered by 
the excluded region from the LEP Higgs searches. The region Ax^ < 4.61 could be accessed 
already with 2 fb~^ in case P3, and 8 fb^^ would give access to the region with Ax^ < 2.30. 
However, even the Ax^ < 4.61 region of the P4 surface would be inaccessible with 8 fb . 

^Three other recent evaluations yield slightly different numbers [66-69], but similar discrepancies with 
the SM prediction. 

^ We note that both CDF and DO have already recorded more than 2.5 fb^^ of integrated luminosity. 
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Figure 3: The same (M^,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) P2, 
(c) P3 and (d) P4 as in Fig. [21 displaying also the expected 95% C.L. exclusion sensitivities 
of searches for H/A t^t" at the Tevatron collider with 2, 4, 8 fb~^ in each of the CDF 
and DO experiments (see text). 
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We note that the CDF Collaboration has recently reported a ~ 2-a excess of candidate 
H/A T^T~ events [9], which would correspond to Ma ~ 160 GeV and tan/5 > 45. As 
discussed in Ref. [30], taking into account all the available experimental constraints, this 
possible excess could be accommodated within the NUHM only for rather different values of 
the parameters from those considered in the benchmark scenarios, namely mi/2 ~ 650 GeV, 
mo ~ 1000 GeV, Aq ~ —1900 GeV, ^ ~ 385 GeV. A likelihood analysis yields values of 

~ 9-10, somewhat higher than the values for the benchmark surfaces. Within the four 
benchmark scenarios here, the precision observables are not in good agreement with low 
and large tan /3, reflecting the fact that the points with Ma ~ 160 GeV and tan /? > 45 lie 
well outside the regions with Ax^ < 4.61 on all of these benchmark surfaces. 

5 LHC Phenomenology 

In this Section we present and compare the sensitivities of various LHC searches for MSSM 
Higgs bosons as functions of Ma and tan/5 in the benchmark surfaces PI, P2, P3 and P4. 
The Higgs bosons can either be produced 'directly' or via cascades, starting with gluino or 
squark production [71]. We focus here on the first possibility, but it should be kept in mind 
that the production via cascades could offer additional channels for the Higgs detection. A 
full evaluation of these channels across the benchmark surfaces must await a more complete 
evaluation of the experimental sensitivities to such decay modes. 

We start the analysis with the light MSSM Higgs boson that behaves like the SM Higgs 
boson for Ma ^ Mz- As a consequence, the region Ma ^ Mz can be covered in all 
benchmark scenarios if a SM Higgs with Mf^ = Mh is accessible at the LHC [11, 12, 14]. In 
Fig. m we display on the WMAP-compatible (M^,tan/?) planes the 5-cr discovery contours 
for pp ^ h —>■ '-f'-f a.t the LHC with 30 fb~^ in the CMS detector [14], where the areas to 
the right of the lines (i.e. for larger Ma) are covered by the pp ^ h ^ '-fy search. This 
channel is particlarly important for a precise mass measurement of the lightest MSSM Higgs 
boson. We show separately the sensitivities for a cut-based analysis (blue solid line) and for 
an "optimized" analysis (black dotted line), see Ref. [14] for details. The cut-based analysis 
should be regarded as a conservative result, while the "optimized" analysis should perhaps be 
regarded as an optimistic expectation [72]. In the cases of surfaces PI and P2, the LHC cut 
analysis for the /i — > 77 search covers all of the Ax^ < 2.30 region and the optimized 

analysis nearly the whole parameter plane. For P3 only parts of the preferred region can 
be covered, while for P4 even with the optimized analysis the best-fit point as well as large 
parts of the Ax^ < 2.30 area remain uncovered. In this region, more luminosity would need 
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Figure 4: The same (M^,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) 
P2, (c) P3 and (d) P4 as in Fig. [2l displaying the expected sensitivities of searches for 
— » /i ^ 77 at the LHC with 30 fb~^ in the CMS detector using a cut analysis or an 
"optimized" analysis (see text) as well as the searches for W^W~ ^ h ^ t^t^ with 60 fb~^ 
in the CMS detector. The parameter regions to the right of the contours are covered at the 
5-cr level. For PI and P2 the W^W~ — > /i — r+r~ channel covers the whole region of the 
(M^,tan/5) plane that is unexcluded by LEP. 
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to be accumulated in order to see a 5-0" signal in the pp ^ h ^ ^y^y channel. 

We turn next to the reaction W'^W'' h t^t^ . On the WMAP-compatible 
(M4,tan/3) planes in Fig. H] we display the 5-cr discovery contours for W^W~ h ^ t^t~ 
at the LHC with 60 fb"^ in the CMS detector [14], where the areas to the right of the lines 
(i.e. for larger M4) are covered by this search. In the cases of surfaces PI and P2, the 
5-cr discovery contours lie within the region already excluded by LEP, so this search cov- 
ers all the unexcluded parts of the surfaces. In the cases of surfaces P3 and P4, however, 
the W^W~ —>■ h —>■ r^T^ discovery contours leave uncovered narrow strips at low Ma for 
tan/3 > 11, 14, respectively. In this part of the parameter space the search for H t^t~ 
should be investigated. In all cases, the 5-0" discovery contours cover the entire Ax^ < 4.61 
regions. However, we note that this channel does not permit a very accurate measurement 
of Mfi, unlike the pp h '-fy channel. 

We now turn to the heavy MSSM Higgs bosons. In Fig. [5] we display in the (MA,tan/3) 
planes the 5-0" discovery contours for bb — > H/A t^t~ at the LHC, where the r's decay 
to jets and electrons or muons (in the BR evaluation for the heavy Higgs bosons possible 
decays to SUSY particles [73-75] have also been taken into account). The analysis is based 
on 60 fb~^ for the final state t^t^ jets [76] and on 30 fb^^ for the r+r^ e+ jet [77] and 
r+r^ fi + iet [78] channels, collected with the CMS detector. As shown in Ref. [79], the 
impact of the supersymmetric parameters other than M4 and tan (3 on the discovery contours 
is relatively small in this channel, and the decays of H/A to SUSY particles [73-75] are in 
general suppressed by large sparticle masses. Only in P4 the decay to the lightest neutralinos 
and charginos is possible over nearly the whole plane (see also Sect. [2]). Including such decays 
in the evaluation of the discovery reach could increase the coverage for heavy Higgs bosons 
somewhat. As a consequence of the relatively small impact of the other SUSY parameters, 
the discovery contours in the four benchmark surfaces are similar to each other and to those 
in the "conventional" benchmark scenarios [79]. The 5-0" discovery contours for the various r 
decay modes are shown separately: they may each be scaled individually for different values 
of the jet (j), /i and electron (e) detection efficiencies, see Ref. [79]. The sensitivities of the 
three different search strategies could in principle be combined, but information required 
for making such a combination is not yet available from the CMS Collaboration. Nor is 
the information available that would be needed to extend the discovery contours to small 
Ma < 200 GeV or to large Ma > 500 to 800 GeV. Nevertheless, we see that the whole 
Ax^ < 2.30 regions of the surfaces PI and P2 would be covered by the LHC H/A t^t~ 
searches, and most of the corresponding regions of the surfaces P3 and P4. Comparing the 
LHC sensitivities shown in Fig. [5] with the Tevatron sensitivities shown in Fig. [3], we see that 
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Figure 5: The same (M^,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) P2, 
(c) P3 and (d) P4 as in Fig. [21 displaying the 5- a discovery contours for H/A t^t~ 
at the LHC with 60 or 30 fb~^ (depending on the r decay channels) and for —>■ r^u 
detection in the CMS detector when Mh± > fnt (see text). 
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the LHC provides access to considerably heavier H/A, up to about 800 GeV, and that the 
covered region extends to lower values of tan/3, reaching tan/5 ~ 10 at low Ma- Comparing 
with Fig. m we see that the H/A t^t^ searches presumably also cover the regions at 
Ma < 150 GeV and tan/3 > 11,14 that were left uncovered in the P3 and P4 surfaces, 
respectively, by the W~^W~ —>■ h ^ r^r^ searches. It would be interesting to verify this by 
means of an extension of the available CMS analysis. 

We also show in Fig. [5] the 5-cr contours for discovery of the via its t^u decay mode 
at the LHC, in the case Mh± > rrit. We see that the coverage is limited in each of the 
scenarios PI, P2, P3 and P4 to Ma < 300 GeV and tan/3 > 30, reaching a small part of 
the Ax^ < 2.30 region of surface P3, only a small part of the Ax^ < 4.61 region of surface 
PI, and not even reaching this region in scenarios P2 and P4. One may also search for 
T^v for lighter Mh± < rrit, but in the cases of surfaces PI and P2 this would be 
useful only in the regions already excluded by LEP, and the accessible regions in surfaces 
P3 and P4 would also be quite limited. 

Another class of possible measurements at the LHC comprises the precise determinations 
of h decay branching ratios [80], and using their ratios to search for deviations from the SM 
predictions for a Higgs boson of the same mass. Such deviations may arise in the MSSM due 
to differences in the tree-level couplings and due to additional (loop) corrections. The most 
sensitive observable is likely to be the ratio of BR(/i r+r^)/BR(/i — ^ VTiy*). We display 
in Fig. [6] the 1-, 2-, 3- and 5- a contours (2- cr in bold) for SUSY induced deviations of this 
ratio of branching ratios from the SM prediction (with Mf^ = Mh). The contours correspond 
to an integrated luminosity at the LHC of 30 or 300 fb~^ [81] (assuming SM decay rates). An 
experimental resolution for BR(/i ^ r+r-)/BR(/i ^ WW*) between 30% (28%) and 45% 
(33%) can be achieved for 30 (300) fb~^. For Mh = 120 GeV the corresponding precision is 
38% (29%). The most promising surfaces are P3 and P4, and we see that over essentially 
all the left lobe of the Ax^ < 4.61 region for P4 a 5-cr discrepancy with the SM should be 
detectable^. On the other hand, only partial coverage of the left lobe of surface P3 would be 
possible, and the sensitivities in the right lobes of P4 and P3 and in the PI and P2 surfaces 
are considerably less promising. Nevertheless, measuring BR(/i r+r^)/BR(/i WW*) 
does offer the prospect of distinguishing between the NUHM and the SM in the low Ma 
regions of surfaces P3 and P4. 



^ It should be kept in mind that the actual experimental precision on the ratio BR(/i ^ r+r )/BR(/i - 
T4^M^*) will be different in this parameter region from the numbers quoted above which assume SM rates. 
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Figure 6: The same (M^,tan/3) planes for the NUHM benchmark surfaces (a) PI, (b) P2, 
(c) P3 and (d) P4 as in Fig. [21 displaying the 1-, 2-, 3- and 5- a contours (2-cr in bold) for 
SUSY-induced deviations on the ratio BR{h r+r-)/BR(/i WW*) at the LHC with 
30 or 300 fb~^ (see text). In the case of surface P2, only 1-a curves are seen in the lower 
part of the figure. The upper curves correspond to 0- a. 
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6 ILC Phenomenology 



In this section we analyze the deviations in the branching ratios of the hghtest MSSM Higgs 
boson to SM fermions and gauge bosons in comparison with a SM Higgs boson of the same 
mass that could be measured at the ILC (see also Ref. [21]). The experimental precisions 
for the branching ratios we analyze are summarized in Tab. [2l 



collider 


channel 


exp. precision [%] 


ILC(500) 


BR{h bb) 


1.5 


ILC(500) 


BR{h T+T-) 


4.5 


ILC(500) 


BR{h WW*) 


3.0 


ILC(IOOO) 


BR{h bb)/BR{h WW*) 


1.5 



Table 2: Experimental precisions at the ILC for various branching ratios of the lightest 
MSSM Higgs boson (assuming SM decay rates) [18,82,83]. The experimental precision in 
the last column corresponds to 1 o" in the plots below. ILC(500,1000) refers to a center-of- 
mass energy of 500, 1000 GeV, respectively. 

We show in Fig. [7] the prospective sensitivity of an ILC measurement of the BR{h —>■ bb) 
in the four (Myi,tan/3) planes. The experimental precision is anticipated to be 1.5%, see 
Tab. [21 We display as solid (blue) lines the contours of the +5, +3, +2, +1, cr deviations 
(with +2 cr in bold) of the MSSM result from the corresponding SM result (for low and 
large tan/? in P2 we also find contours for —2, —1 a, with — 2cr in bold). The separations 
between the contours indicate how sensitively the SUSY results depend on variations of Ma 
and tan (3. Also shown in Fig. [7] via dashed (green) lines is the sensitivity to SUSY effects 
of the ILC measurement of the ratio of branching ratios BR{h — » bb)/BR{h WW*) (for 
low Ma and large tan/? in P2 we also find contours for —5, —3, —2, —la). The precision 
measurement of the ratio BR{h —>■ bb)/BR{h WW*) clearly provides a much higher 
sensitivity to SUSY effects than the measurement of BR{h — > bb) alone (see also Ref. [20]). 

For the ILC measurement of the BR{h bb), in the cases of PI and P2 we see that 
the prospective sensitivities are less than 3 a throughout almost all the regions with Ax^ < 
4.61. The situations are different, however, for the planes P3 and P4. In each case, the 
cosmologically-favoured region is divided into separate lobes at low and high Ma- In the 
P3 case, the measurement of BR{h bb) would be sufficient to establish a SUSY effect 
with more than five a throughout most of the I0W-M4 lobe, and all of it in the P4 case. A 
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precision measurement of BR(/i 66)/BR(/i WW*) yields a significant improvement for 
all benchmark surfaces. We see that, in case PI, the sensitivity already exceeds 5 a in much 
of the region with Ax^ < 2.30, and the fraction of this region covered at the 5-cr level is 
even larger in the case P2. Even more encouragingly, in the case P3 the sensitivity exceeds 
5 a throughout the Ax^ < 2.30 region, and in the case P4 it exceeds 5 a by a substantial 
amount throughout the Ax^ < 4.61 region. 

Next, we show in Fig. [8]the prospective sensitivity of an ILC measurement of the BR(/i 
r+r^) in the four (Myi,tan/5) planes, using solid (red) contours. In the cases of PI and P2, 
we again see that the prospective sensitivities are less than 3 a throughout almost all the 
regions with A^^ < 4.61. In the cases of planes P3 and P4, the sensitivities are greater, but 
less than the corresponding sensitivities to the BR(/i —>■ bb) shown previously in Fig. [71 Of 
all the single ILC measurements, the one with the greatest sensitivity to SUSY effects is that 
of the BR(/i —>■ WW*), which is also shown in Fig. [8] using dashed (black) lines. In the cases 
PI and P2, we see that the sensitivity may rise above 5 a already within the Ax^ < 4.61 
region. In the case of P3, the sensitivity is well above 5 cr throughout the low-M^ region. 
In the case of P4, a 5-cr significance is exceeded already in much of the high-M^ lobe, where 
the sensitivity never falls as low as 3 cr in the favored region. 

We have not made a complete study of the combined sensitivity of the ILC measurements 
to the benchmark surfaces, but it is clear from this brief survey that the ILC measurements 
would in general provide interesting tests of the MSSM at the loop level. In the absence of 
detailed studies, we expect that CLIC measurements would have similar sensitivities, since 
h production would be more copious at the higher CLIC energies, and the CLIC luminosity 
at lower energies could be similar to that of the ILC [23]. In addition to the precision 
measurements described here, the ILC and CLIC would be able to produce directly associated 
H + A pairs above the kinematic threshold. 

7 B Physics 

We display in Fig. [9] the results for three BPO BR{b 57), BR(S, ^ /i+yU"), BR(S„ ^ 
TUr), in the four benchmark {Ma, tan (3) planes. 

The prediction of /^"""/i" is based on Ref. [47,84]. The solid (beige) line indicates 

BR(i?s yU'^'/U") = 10"'', corresponding roughly to the current upper bound from CDF [85] 
and DO [86]. The latest bound reported by CDF has recently been lowered to 5.8 x 10~^ [87]. 
The dashed (beige) line indicates a BR of 2 x 10"^. In Fig. [9] we see that the current upper 
limit on jj,'^ ji~ already excludes regions of the planes at small Ma and large tan/?. 
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Figure 7: The same (M^,tan/3) planes for the NUHM benchmark surfaces (a) PI, (b) P2, 
(c) P3 and (d) P4 as in Fig. displaying 5, 3, 2, 1, O-cx sensitivity contours (2-cr in bold) for 
SUSY effects on BR{h bb) (solid blue lines) and BR{h bb)/BR{h WW*) (dashed 
green lines) at the ILC (see text). Note that for surface P2 for low Ma and large tan/3 
also —2,—1-a are shown for BR(/i — > bb), and — 5, — 3, — 2, — 1-cr are shown for BR(/i — >• 
bb)/BR{h WW*). 
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Figure 8: The same (MA,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) P2, 
(c) P3 and (d) P4 as in Fig. [21 displaying 5, 3, 2, l-cr sensitivity contours for SUSY effects 
on the BR(/i t^t^) at the ILC (solid red lines). Also shown are the —5, —3, —2, —l-cr 
sensitivity contours for the SUSY effects on BIi{h — > lyVT*) at the ILC (dashed black lines). 
Note that for surface P2 =b2, ±1 and 0-a are shown for BR(/i —>■ r+r^), and ±5, ±3, ±2, ±1 
and 0-a are shown for BR(/i — > iyiy*). 



21 



starting to cut into the region with A^^ < 4.61. The prospective sensitivities would extend 
as far as the best-fit points. 

For b —* s'y our numerical results have been derived with the BR(6 57) evaluation 
provided in Refs. [88], incorporating also the latest SM corrections provided in Ref. [89]. 
The results in Fig. [9] are shown as the two blue lines indicating BR(6 — > 57) of 4 x 10"'* 
(solid) and 3 x 10~^ (dashed). These have to be compared to the experimentally preferred 
value of BR(6 ^ 57) = (3.55 ± 0.241°;°^ ± 0.03) x 10"^ [90]. The best-fit point together with 
large parts of the preferred regions lie between the two lines, i.e., large parts of the four 
benchmark planes are in good agreement with the current experimental value. 

Our results for BK{Bu tVt) are based on Ref. [91]. In the four benchmark scenarios 
of Fig. [9] the results are shown in form of the NUHM result divided by the SM prediction 
as black lines. The solid (dashed) lines correspond to a ratio of 0.9 (0.7), where the current 
central value is 0.93 ±0.41 [92,93]. It can be seen that the best fit value as well as large parts 
of the preferred parts of the benchmark planes predict a value somewhat lower than the 
current experimental result. However, with the current precision no firm conclusion can be 
drawn. 

8 Direct Detection of Supersymmetric Dark Matter 

In Fig. [To] we show how the direct detection of the LSP via spin-independent scattering on 
nuclei probes the four (M^,tan/3) planes. We focus here on the bound from the XENONIO 
experiment that was recently published by the XENON collaboration [27], which improves 
on the previous CDMS results [26]. We note that the XENONIO experiment has seen some 
potential signal events which are, however, interpreted as background. 

The constraint imposed by the limits from direct detection experiments is sensitive to two 
theoretical uncertainties that are independent of the specific dark matter model. One is the 
local density of cold dark matter, which is normally estimated to be Pcdm = 0.3 GeV/cm'^, 
although smaller values may be consistent with some models of the Galaxy. The other im- 
portant uncertainty is that in the nucleonic matrix element of the local operator responsible 
for the spin-independent scattering amplitude. This is related, in particular, to the so-called 
a term, S^jy? that may be derived from measurements of low-energy vr-nucleon scattering. 

The solid lines in Fig. [10] correspond to the XENONIO bound obtained assuming Pcdm = 
0.3 GeV/cm^ and using St^at = 45 MeV as input, corresponding to a relative strange-quark 
density y = 2{N\ss\N) / {N\{uu + dd)\N) = 0.2 [94]. These assumptions are realistic, though 
there is a large uncertainty in the strangeness contribution which may lead to larger rates if 
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Figure 9: The same (M^,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) 
P2, (c) P3 and (d) P4 as in Fig. [21 displaying the expected sensitivities of the B physics 
observables Bg —>■ fi'^fi~, b —>■ s'-f and B^ —>■ rv. The various lines indicate: BR(i?5 

= 10-^(2 X 10"^) as solid (dashed), BR(6 ^ 57) = 4(3) x 10"^ as solid (dashed), 
BR(i?„ t"z/^)mssm/sm = 0.9(0.7) as sohd (dashed). 
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Figure 10: The same (Myi,tan/5) planes for the NUHM benchmark surfaces (a) PI, (b) 
P2, (c) P3 and (d) P4 as in Fig. [2], displaying the expected sensitivities of present and 
prospective future direct searches for the scattering of dark matter particles (see text). 
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S^AT is larger or significantly lower rates if the strangeness contribution to the proton mass 
is small. The dashed lines show the bounds that one would obtain from the XENONIO 
experiment assuming the same value of Pcdm, but with E^at = 36 MeV corresponding to 
y = 0, and therefore representing more conservative assumptions. Finally, as an example of 
the possible sensitivity of future experiments, the dotted lines show the contours one would 
obtain for a spin-independent cross section of 10^*^ pb, assuming the same value of pcdm and 
S^AT = 45 MeV as input. 

We see from Fig. [10] that the surfaces PI and P2 are not probed by the current limits 
from the XENONIO experiment. Only the possible future sensitivity at 10^® pb starts to cut 
into the Ax^ < 4.61 region. For these planes, accelerator searches are clearly more powerful. 
The situation is different for the planes P3 and P4, due to the relatively low values of mi/2 
across these planes. We recall that, for planes PI and P2, mi/2 scales with and the 
sparticle spectrum is typically heavier at large Ma than at the corresponding points in planes 
P3 and P4. As a result, the spin-independent Xi ~ P elastic cross section is suppressed for 
planes PI and P2. On the other hand, we see that the current XENONIO bound probes 
large parts of the Ax^ < 2.30 areas of P3 and P4 planes, if one uses the moderate values 
of Sttat = 45 MeV and the strange-quark content. Indeed, in the case of the P4 surface, the 
current XENONIO bound would even cover the best-fit point for this value of S^^r and the 
default value for the local density of cold dark matter. The more conservative analysis, on 
the other hand, is sensitive only to smaller Ma values, and probes only a much smaller part 
of the regions preferred by the analysis. Finally, we note that a future sensitivity to a 
cross section of 10"^ pb would cover the entire PS and P4 surfaces. 

9 Conclusions 

The value of benchmark studies is that they allow one to understand better the range of 
possibilities opened up by supersymmetry. It is therefore desirable that benchmarks be 
chosen in such a way as to respect, as far as possible, the definitive experimental constraints, 
and also that they be susceptible to systematic study. We have demonstrated in this paper 
how NUHM benchmark surfaces chosen so that the relic cold dark matter density falls 
within or below the range favoured by WMAP and other experiments may be used to probe 
supersymmetric phenomenology. Our approach based on the NUHM scenario significantly 
differs from previous proposals of benchmark scenarios for the MSSM Higgs sector that 
were entirely formulated in terms of low-scale parameters and that were not suitable for 
a phenomenologically acceptable prediction of the cold dark matter density. The analysis 
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of our proposed benchmark surfaces is facilitated by developments in the FeynHiggs code 
that are described in the Appendix. These will enable the interested reader to explore the 
prospects for her/his favourite experimental probe of supersymmetry in these benchmark 
surfaces. 

We have displayed the constraints currently imposed in the new benchmark surfaces 
by electroweak precision observables, and explored the prospects for Higgs searches at the 
Tevatron coUider, the LHC and the ILC, and we have also explored indirect effects in B 
physics and in dark matter detection. Whereas the Tevatron collider will be able only to 
nibble at corners of these NUHM benchmark surfaces, experiments at the LHC will be able 
to cover them entirely, and the ILC will have good prospects for precision measurements. 
There are good prospects for B experiments in parts of the benchmark surfaces, and direct 
dark matter may be detectable in some cases. 

It should of course be noted that benchmark studies may soon be rendered obsolete - 
namely by the discovery of supersymmetry. 

As we were completing this paper, we heard the sad news of the passing away of Julius 
Wess, one of the discoverers and founding fathers of supersymmetry. Julius did so much to 
develop our understanding of supersymmetry, to awaken our appreciation of its beauty, and 
to convince us of its importance for physics. Humbly and respectfully, we dedicate this paper 
to his memory. 
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A Evaluation of Benchmark Surfaces with FeynHiggs 



The new benchmark surfaces have been implemented into the code FeynHiggs [54-57]. In 
this way, any user may apply them to perform phenomenological analyses. 

From the mathematical point of view, the NUHM/CDM constraints introduce non-trivial 
relations between input parameters, which thus cannot be scanned naively by independent 
loops. To solve this in a generic way, FeynHiggs 2.6 allows the user to interpolate the inputs 
from a Parameter Table into which arbitrary relations can be encoded. The tables containt- 



ing the four benchmark surfaces can be downloaded from |http : / /www . f eynhig gs . de . To 
implement the new format of a Parameter Table, significant internal rearrangements were 
necessary from which the concept of a FeynHiggs Record evolved. 

A Record is a new data type which captures the entire content of a parameter file in 
the native format of FeynHiggs. In this respect it is akin to the SUSY Les Houches Accord 
Record [95], but also encodes information about parameter loops and has 'inheritance rules' 
for default values. Using the routines to manipulate a Record, the programmer can, among 
other things, process FeynHiggs parameter files independently of the front end. 

In addition to containing loops over parameters, a Record can be associated with a 
Parameter Table in such a way that values not explicitly given in the parameter file are 
interpolated from the table (as it can be done for the four benchmark scenarios). 

The FeynHiggs Record is conceptually a superstructure 'on top' of the conventional part 
of FeynHiggs. This means that a Record can be manipulated without any influence on the 
computation of Higgs observables at first. Only when the FHSetRecord subroutine is invoked 
are its current values set as the inputs for the computation. So in principle, the FeynHiggs 
Record can be used without doing any computation of Higgs observables at all. 

Technically, a Record is a two-dimensional real array of the form 



rec(?|,j^) 


iVar 


iLower 


iUpper 


iStep 


iTB 


L 


U 


U 


U 


iMAO 


L 


U 


U 


U 



The column index i specifies the parameter. The indices are labelled as in the parameter 
file, but prefixed with an i (see Table |3]). 

The row index j enumerates the variables that constitute the loop over a parameter, 
i.e. the current, lower, and upper value and the step size. The loop inferred through 
these parameters has the form 
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do rec(i,iVar) = rec(i,iLower) , rec(i,iUpper) , rec(i,iStep) 
enddo 

• U entries indicate fields filled in by the user. If no loop is desired over a particular 
parameter, the fields rec(z,iUpper) and rec(z,iStep) can be omitted. On top of 
that there are also 'inheritance rules' (given in Table [3]), stating for example that M3SL 
defaults to MSusy if not given explicitly. 

• L entries indicate fields replaced by the FHLoopRecord routine while working off the 
loops over parameter space, i.e. these fields are updated automatically according to the 
current point in the loop. For example, if the Record contains 

rec(iTB, iLower) = 10 
rec(iTB, iUpper) = 50 
rec(iTB,iStep) = 10 

the first call to FHLoopRecord will set rec(iTB,iVar) to 10, the next to 20, etc. 

A.l Fortran Use 
A. 1.1 Declaration 

Every subroutine or function which uses a Record must first include the definitions: 
#include "FHRecord.h" 

Records can then be declared with the preprocessor macro RecordDecl, which hides the 
declaration details. For example, 

RecordDecl (rec) 

declares the Record rec (i ,j) . When declaring several records, each needs its own RecordDecl 
statement, i.e. RecordDecl (reel , rec2, . . .) is not permissible. 

A. 1.2 Initializing a Record 

A FeynHiggs Record has to be brought into a defined state before its first use, either by 
clearing it with 

call FHClearRecord(rec) 
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or by reading it from a file, whicfi similarly overwrites any previous content: 

call FHReadRecord (error , rec, "file") 
if( error .ne. ) stop 

where file is the name of a parameter file in FeynHiggs' native format. 
Fields can be set or read out using ordinary Fortran array access, e.g. 

rec(iTB, iLower) = 10 
or print *, "At = ", rec(Re(iAt) , iVar) , rec(lni(iAt) ,iVar) 

The 'current value' field (iVar) should not be set explicitly, as it is updated automatically 
by FHLoopRecord. 

A. 1.3 Looping over a Record / Setting the FeynHiggs input 

The loops over parameters contained in a Record are worked off through calls to FHLoopRecord, 
which update the Record's 'current value' fields (iVar). FHLoopRecord is thus usually in- 
voked in the context of a looping construct, such as 

call FHLoopRecord (error , rec) 
do while ( error .eq. ) 

call FHLoopRecord (error , rec) 
enddo 

The subroutine FHSetRecord can be used to set the 'current value' fields (iVar) as input 
parameters for FeynHiggs. This works effectively as a combination of FHSetPara, FHSetCKM, 
and FHSetNMFV, except that the parameters are taken from the Record. In a typical appli- 
cation the above loop would be extended to 

call FHLoopRecord(error , rec) 
do while ( error .eq. ) 

call FHSetRecord (error , rec, IDO) 

if( error .ne. ) stop 

call FHHiggsCorr (error , MHiggs, SAeff, UHiggs, ZHiggs) 
if( error .ne. ) stop 

call FHLoopRecord (error , rec) 
enddo 
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The third argument in FHSetRecord is the same scale factor which appears in FHSetPara 
and which determines the renormahzation scale as a multiple of the top mass. 

A. 1.4 Associating a Record with a Table 

The FeynHiggs Record allows one to interpolate parameters from a data table. The table is 
interpolated in two user-selectable variables which can be chosen identical if interpolation in 
only one variable is desired. 

The table first needs to be loaded into internal storage. At the moment FeynHiggs has 
a static allocation for one table of at most 2400 hues. This allows the complete implemen- 
tation in Fortran and seems sufficient for all present applications. The table's format is 
rather straightforward: the first line contains the column names (same identifiers as in the 
FeynHiggs input file), followed by the data rows. All items are separated by whitespace. 

Loading the table can either be done through the input file and is thus automatically 
performed in FHReadRecord. To this end one has to add a line 

table file varl var2 

to the parameter file. For example, "table mytable TB MAO" reads the file mytable into 
memory and sets TB and MAO as input variables for the interpolation. The table must 
obviously contain columns for the input variables. 

It is also possible to integrate the table file into the parameter file. The table statement 
then takes the form 

table - varl var2 

and must be the last statement in the parameter file, followed immediately by the table data. 
Alternately, the table is loaded by 

call FHLoadTable (error, "file", 5) 
if( error .ne. ) stop 

The table is read from file, unless that equals in which case the table is read from the 
Fortran unit given in the third argument (unit 5 is Fortran's equivalent of stdin and hence 
a good default argument here). 

The table is associated with the record through 

call FHTableRecord (error , rec, varl, var2) 
if( error .ne. ) stop 
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where varl and var2 are the indices of the input variables, e.g. iTB and iMAO. To translate 
parameter names (strings) into indices, one can use the FHRecordlndex subroutine, as in: 

call FHRecordlndex (index, name) 

A. 2 Mathematica Use 

Using FeynHiggs Records in Mathematica is for the larger part very similar to doing so in 
Fortran. The main difference is that one does not have to declare a Record. Rather, both 
initialization routines 'create' the Record: 

rec = FHClearRecord [] 
or rec = FHReadRecord["f ile"] 

The Record is represented as an FHRecord object in Mathematica. Access to fields is very 
similar to the Fortran case, e.g. 

rec[[iTB,iLower]] = 10 
or Print ["At = rec [ [Re [iAt] , iVar] ] , rec [ [Iin[iAt] , iVar] ] ] 

So is the use of FHLoopRecord, except that the updated Record is returned, rather than 
modified in situ. In other words, FHLoopRecord returns an FHRecord as long as the loop 
continues. The loop would thus look hke 

While [ Head [rec = FHLoopRecord [rec]] === FHRecord, 
] 

The other routines are used straightforwardly, for example: 

FHSetRecord[rec, 1] 

FHLoadTable["file"] 

rec = FHTableRecord[rec, varl, var2] 

index = FHRecordlndex [name] 

A. 3 Examples 

A. 3.1 Command-line Mode with Peirameter Table 

In the simplest case, a Parameter Table can be processed through an input file with a table 
statement: 
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MAO 203 
TB 5.7 

table file.dat MAO TB 

The Parameter Table is read from file.dat in a format like 



MT 




MSusy 


MAO 


TB 


At 


MUE 


171 


.4 


500 


200 


5 


1000 


761 


171, 


.4 


500 


210 


5 


1000 


753 


171, 


.4 


500 


200 


6 


1000 


742 


171, 


.4 


500 


210 


6 


1000 


735 



Alternately, the Table can be integrated into the parameter file, as in 

MAO 203 
TB 5.7 
table - MAO TB 
MT MSusy MAO TB 
171.4 500 200 5 
171.4 500 210 5 

171.4 500 200 6 
171.4 500 210 6 

This minimal setup assumes that all parameters are contained in the table. More generally, 
the ones not contained in the table have to be given in the parameter file. The interpolation 
for the parameters given (here MAO and TB) is performed automatically by FeynHiggs. 

A. 3. 2 Using a Record with Table in Fortran 

In Fortran, the same example might be coded as 

program record_test 
implicit none 

#include "FHRecord.h" 



At MUE . . . 

1000 761 

1000 753 

1000 742 

1000 735 
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RecordDecl (rec) 
integer error 

double precision MHiggs(4) 

double complex SAeff, UHiggs(3,3), ZHiggs(3,3) 

call FHClearRecord(rec) 
rec(iMAO,iLower) = 203 
rec(iTB,iLower) = 5.7 

call FHLoadTable (error, "file.dat", 5) 
if( error .ne. ) stop 

call FHTableRecord (error, rec, iTB, iMAO) 
if( error .ne. ) stop 

call FHSetFlags(4, 0, 0, 3, 0, 2, 1, 1, 3) 

call FHLoopRecord (error , rec) 
do while ( error .eq. ) 

call FHSetRecord (error , rec, IDO) 

if( error .ne. ) stop 

call FHHiggsCorr (error, MHiggs, SAeff, UHiggs, ZHiggs) 
if( error .ne. ) stop 

print *, "TB, Mhl = ", rec(iTB, iVar) , MHiggs(l) 

call FHLoopRecord (error , rec) 
enddo 
end 

A. 3. 3 Using a Record with Table in Mathematica 

In Mathematica, the structure and syntax is very similar to Fortran (mainly round brackets 
have to be converted into square ones): 
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Install ["MFeynHiggs"] 

rec = FHClearRecord[] 

rec[[iMAO,iLower]] = 203; 
rec [ [iTB , iLower] ] =5.7 

FHLoadTable [ " f ile . dat " ] 

rec = FHTableRecordErec, iTB, iMAO] 

FHSetFlags[4, 0, 0, 3, 0, 2, 1, 1, 3] 

While [ Head [rec = FHLoopRecord[rec]] === FHRecord, 
FHSetRecord[rec, 1]; 
res = FHHiggsCorr [] ; 

Print["TB, Mhl = rec [ [iTB, iVar] ] , (MHiggs /. res)[[l]] ]; 

] 
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Table 3: The parameter index names of a FeynHiggs Record. Indices of real parameters are 
listed in the left, of complex ones in the right column. Complex quantities, e.g. Af, can be 
accessed either through Re (iAt) and Im(iAt), or Abs(iAt) and Arg(iAt), with iAt alone 

as a synonym for Re (iAt). In cases where both Re/Im and Abs/Arg are given, the latter 
take precendence. Please consult the FeynHiggs(l) manual page for more details. 



Inck'x iiamc^ 




Dcfaull value 


lAlfasMZ 


«s(M|) 


-1 


IMC 


nic 


-1 


IMT 


mt 




IMB 


mb(on-shell) 


-1 


IMW 


Mw 


-1 


iMZ 


A /f 

Mz 


— 1 


TB 


tan B 




MAO 






MHp 


Mh+ 




IMSusy 


MsvSY 




IMSSL 


M| 


IMSusy 


iM2SL 


M| 


iM3SL 


IMISL 


Mi 


iM2SL 


iM3SE 


M| 


IMSusy 


iM2SE 


M| 


iM3SE 


IMISE 


M| 


iM2SE 


iM3SQ 


IMSusy 


iM2SQ 


M? 
Q 


iM3S0 


iMloLj 


!\/f 1 
f 




iM3SU 




IMSusy 


iM2SU 


M? 
u 


iM3SU 


IMISU 


Ml 

4 


iM2SU 


iMSSD 


IMSusy 


iM2SD 


Ml 

..^i 


iM3SD 


IMISD 


iM2SD 


IQtau 


Qr 





iQt 


Qt 





IQb 


Qb 





iCKMthetal2 


012 


-1 


iCKMtheta23 


023 


-1 


iCKMthetal3 


013 


-1 


iCKMdeltai3 


Sl3 


-1 



Index name 


r^u'amel er 


Defaiill value 


IMl 


Ml 





iM2 


M2 




iM3 


Ms 




iAt 


At 




iAc 


Ac 


iAt 


iAu 


Au 


iAc 


H AH 




-i A-t- 


iAs 


As 


iAb 


iAd 


Ad 


iAs 


iAtau 


Ar 


iAb 


iAmu 




iAtau 


iAe 


Ae 


iAmu 


ideltaLLuc 







ideltaLRuc 


^LR 
^uc 





ideltaRLuc 


^RL 





ideltaRRuc 


^RR 

uc 





ideltaLLct 


.LL 
^ct 





ideltaLRct 


"ct 





ideltaRLct 


xRL 


(J 


ideltaRRct 







ideltaLLut 


Silt 





ideltaLRut 







ideltaRLut 


e 





ideltaRRut 







ideltaLLds 


^ds 





ideltaLRds 







ideltaRLds 







ideltaRRds 


xRR 
"ds 





ideltaLLsb 


"sb 





ideltaLRsb 


"sb 





ideltaRLsb 







ideltaRRsb 


ST 





ideltaLLdb 


"db 





ideltaLRdb 


:LR 
"db 





ideltaRLdb 


ST 





ideltaRRdb 


ST 
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